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ABSTRACT 

We present multi- frequency radio-continuum and Hi observations of NGC 4438, the 
highly-disturbed, active galaxy in the Virgo cluster, with the Very Large Array (VLA) 
and the Giant Metrewave Radio Telescope (GMRT). High-resolution observations of 
the central 1 kpc with the VLA at 4860 and 8460 MHz show the presence of an 
inverted-spectrum radio nucleus located between the highly asymmetric lobes of radio 
emission. This demonstrates that these lobes which are seen in radio continuum. Ha 
and x-ray wavelengths and are located at ~230 and 730 pc from the nucleus arise due 
to an active galactic nucleus (AGN) rather than a compact nuclear starburst. The low- 
frequency radio continuum observations made with the GMRT detect the extended 
emission on the western side of the galaxy whose spectral index is flatter at higher 
frequencies and suggests that it is a mixture of thermal and non-thermal emission. 

The Hi observations show an elongated structure which is displaced by ~4.1kpc 
on the western side of NGC 4438 and has a size of ~9.8 kpc and a mass of 1.8x lO^M©. 
The velocity field suggests systematic rotation, consistent with earlier observations. 
These observations also detect Hi emission from the disk of the galaxy with a mass of 
1.2x lO^M©. We detect a faint Hi-tail towards the north of NGC 4438 close to a stellar 
tail seen earlier in deep optical observations. This Hi-tail has a total extent of ~50 kpc 
and a mass of 1.4x 10® M© if it is at the distance of NGC 4438. The velocity of the Hi 
tail is ~— 10 km s^^ similar to that of Hi emission from IC 3355, but the possibility 
that the tail could be foreground Galactic emission cannot be ruled out. We discuss 
the different structures in the light of different models which have been suggested for 
this disturbed galaxy, namely ram pressure stripping, tidal and ISM-ISM interactions. 

Key words: galaxies: individual: NGC 4438 - galaxies: individual: IC 3355 - galaxies: 
nuclei - galaxies: interactions - galaxies: ISM - galaxies: kinematics and dynamics 



1 INTRODUCTION 

The galaxy, NGC 4438 (VV188, Arp 120), which is located 
only about 1° (~300 kpc) from the centre of the Virgo clus- 
ter, is highly inclined (80°) and has a very disturbed disk 
(Fig. 1) with various components of the interstellar medium 
(ISM) being visible on the western side of the disk. The basic 
properties of the galaxy are summarised in Table 1. High- 
resolution x-ray observations with the Chandra telescope 
show emission from a ~700 pc nuclear region, a ~2.3 kpc 
spherical bulge and a network of filaments extending 4—10 
kpc to the west and south-west of the galaxy (Machacek, 
Jones & Forman 2004). These regions are well correlated 
with similar features seen in Ha (Kenney & Yale 2002; 
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Chemin et al. 2005). Diffuse extended radio emission at 1.4 
GHz has been seen extending up to ~10 kpc on the western 
side of the galaxy. Some have suggested that these features 
on the western side might be due to interactions with the 
intracluster medium (ICM) of the Virgo cluster since the 
relative velocity of NGC 4438 is ~1000 km s"^ A tidal en- 
counter with the SBO galaxy NGC 4435 at a projected dis- 
tance of ~4.3 arcmin (~20 kpc) is also likely to have affected 
NGC 4438. Kenney et al. (1995) have suggested that most 
of the features of the disturbed ISM are likely due to a high- 
velocity ISM-ISM coUision between NGC 4438 and 4435. 
Numerical simulations suggest that the centres of these two 
galaxies have passed within 5—10 kpc of each other ^10* yr 
ago (VoUmer et al. 2005). Radio continuum observations of 
the nuclear region by Hummel & Saikia (1991) have shown 
two well-defined lobes of radio emission in this galaxy which 
is unusual, though not unique, for a spiral galaxy. These 
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Table 1. Basic data on NGC 4438." 



RA (J2000)'' 


DEC (J2000)'' 


Type= 


a 


X h'^ 




if 




D'' q-factor* 


HP 


(h m s) 






(' 


X ') 


O 


O 


(km s-l) 


(Mpc) 


deficiency 


12 27 45.67 


+13 00 31.5 


SA(s)0/a pec 


8.5 


X 3.2 


29 


80 


71±3 


17 1.83 


>1.0 



a Taken from the NASA Extragalactic Database (NED) unless stated otherwise. 

b Position of the radio nucleus from our high-resolution, VLA A-axray image at 8460 MHz. 

c Morpliological type. 

d Optical major and minor axes. 

e Position angle (PA) of major axis from Kenney et al. (1995). 
/ Inclination angle from Kenney et al. (1995). 
g Heliocentric systemic velocity. 

h Distance of NGC 4438 taken drom VoUmer et al. (2005). For this distance 1"=82 pc. 
i q-fax;tor: logarithmic ratio of FIR to radio luminosity (Reddy & Yun 2004). 
j Hi deficiency parameter from Cayatte et al. (1994). 
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Figure 1. The DSS blue-band image of NGC 4438 and its com- 
panion, NGC 4435. The -|- sign marks the position of the radio 
nucleus discussed in this paper. 



lobes arc highly asymmetric with the north-western one be- 
ing much brighter than the south-eastern one. x-ray and 
Ha-|-[Nii] emission are also seen from these two lobes which 
exhibit a bubble- like structure and similar asymmetry in lo- 
cation, brightness and size at these wavelengths (Kenney & 
Yale 2002; Machacek et al. 2004). 

In this paper, we first present our results of radio contin- 
uum observations of the nuclear region and extended emis- 
sion on the western side of the galaxy (Section 3). Wo then 
present the results of H i observations with the GMRT and 
the VLA D-array (Section 4). This is followed by a brief 
discussion and a summary of the results (Section 5). 



2 OBSERVATIONS AND DATA ANALYSIS 

For the continuum data, an observing log for both the 
GMRT and VLA observations as well as some of the ob- 
served parameters of the contirmum images are presented in 
Table 2 which is arranged as follows. Column 1: name of the 
telescope where we also list the configuration for the VLA 
observations. In addition to our own data wo have also anal- 
ysed many sets of archival VLA data on this galaxy. Column 
2: the frequency of the observations; columns 3 and 4: dates 
of the observations and the time, t, spent on the source in 
minutes; columns 5 and 6: the phase calibrator used and its 
flux density estimated from the observations; columns 7, 8 
and 9: major and minor a^xes of the restoring beam in arc- 
sec and its position angle (PA) in degrees; column 10: the 
rms noise in the continuum image in units of mjy/beam; 
columns 11 and 12: the peak brightness and total flux den- 
sity in units of mJy/beam and mjy for the continuum im- 
ages. The total error in the flux density is approximately 5% 
for frequencies higher than 1 GHz and ~10% at 617 MHz. 

The radio continuum observations with the GMRT as 
well as with the VLA were made in the standard fashion 
with each source observation interspersed with observations 
of the phase calibrator. The primary flux density calibra- 
tor was 3C286 or 3C147 whose flux densities were estimated 
on the Baars et al. (1977) scale using the 1999.2 VLA val- 
ues. The bandwidth of the continuum observations with the 
GMRT at 617 MHz was 16 MHz, while for the L-band obser- 
vations it was 8 MHz. The bandwidth for all the VLA contin- 
uum observations was 50 MHz. The data analysis was done 
using the Astronomical Image Processing System (AIPS) of 
the National Radio Astronomy Observatory (NRAO). Since 
GMRT data are acquired in the spectral-line mode with 128 
spectral channels, gain and bandpass solutions were applied 
to each channel before combining them. Self-calibration was 
applied to all the GMRT and VLA data sets. The low- 
resolution 617-MHz map was made with the 3-D cleaning 
option in CLEAN using 16 fields. The self calibrated gains 
were applied to this data set correcting only for phase and 
rejecting all the failed solutions. For all the VLA data sets 
two cycles of phase and one cycle of amplitude self calibrar 
tion were applied except for the VLA A-axray 8460-MHz 
data where only two cycles of phase self calibration were 
applied. 

The information for the Hi observations are presented 
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Figure 2. Left panel: The VLA A-array image at 1490 MHz with an angular resolution of 1.5 arcsec. Contours: 0.037x(— 4, —2.82, 
2.82, 4, 5.65, 8 . . . ) mjy/beain. Middle panel: VLA B-array image at 4860 MHz with an angular resolution of 1.5 arcsec. Contours: 
0.026x(-4, -2.82, 2.82, 4, 5.65, 8 . . . ) mjy/beam. Right panel: Spectral index image between 1490 and 4860 MHz. The contour is at 
—0.8. The + sign marks the position of the radio nucleus in all the images. 
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Figure 3. Maps of the minimum energy parameters of the mini double-lobed radio source, made from Fig. 2, assuming a heavy particle 
to electron energy ratio of k=40. In each case, black denotes the highest values and the star marks the position of the radio nucleus, (a) 
The cosmic ray energy density, ucr. Contours are 200, 350, 750, and 1200 eV cm^"^ and the peak value is 1.58 X 10'^ cV cm^'^ . (b) 
Diffusion length, Lq. Contours arc 100, 150, 200 and 300 pc and the peak is 347 pc. (c) Magnetic field strength, B. Contours arc 50, 65, 
85, and 120 /xG and the peak is 136 /xG. (d) Lifetime of the particles, t. Contours are 1, 2, 2.75, and 4 Myr and the peak is 6.59 Myr. 
(e) Cosmic ray power, P. Contours are 0.33, 0.75, 1.5, and 4 x 10^® ergs and the peak is 5.57 x 10^® ergs s~^. 



© RAS, MNRAS 000, ??-?? 



4 Ananda Hota, D.J. Saikia and Judith A. Irwin 

Table 2. Observation log and observed parameters of the continuum images. 



Telescope 


Frcq. 


Obs. 


t 


Phase 


Seal. 


Beam size 




mis 


Spk 


Stot. 




MHz 


date 


min 


Cahb. 


Jy 


maj. 


min. 


PA 


in,Jy 


m,Jy 


mjy 














II 


// 










(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


GMRT 


617 


2002Janll 


300 


1120+057 


3.57 


46 


23 


149 


4 


151 


311 


VLA-A 


1490 


1986May24 


236 


1252+119 


0.98 


1.50 


1.50 




0.04 


19 


87 


VLA-B 


4860 


1986Jul20 


177 


1252+119 


0.62 


1.50 


1.50 




0.03 


7.3 


33 


VLA-A 


4860 


1986May24 


44 


1252+119 


0.62 


0.41 


0.36 


114 


0.04 


0.5 


12 


VLA-A 


8460 


1998Mayl3 


23 


1236+077 


0.71 


0.21 


0.19 


176 


0.03 


0.4 


7.8 


VLA-AB 


14939 


1991Dec05 


13 


1252+119 


0.52 


0.57 


0.43 


127 


0.19 


1.0 


8.9 




Figure 4. Left panel: The VLA A-array image at 4860 MHz with an angular resolution of ~0.38 arcsec. Contours: 0.042x(— 4, —2.82, 
2.82, 4, 5.65, 8 . . . ) mjy/beam. Right panel: The VLA A-array image at 8460 MHz with an angular resolution of ~0.19 arcsec showing 
only the nucleus and the western shell. Contours: 0.025x(— 4, —2.82, 2.82, 4, 5.65, 8 . . . ) mJy/beam. 



later in Section 4 (Table 5) where, in addition to the values 
described previously, we have also listed the spectral reso- 
lution in units of km s"'^ (column 6), the rms noise in the 
channel maps in units of mjy/bcam (column 10), the rms 
noise in the spectrum in units of mJy (column 11) and the Hi 
flux density in units of Jy km (column 12). The last row 
corresponds to the VLA-D array data to which a broadscale 
taper has been appHod (see Sect. 4.2). 

The analysis of the Hi observations was also done in 
the standard way. 3C286 was the primary flux density and 
bandpass calibrator. The total bandwidth for the GMRT 
H I observations was 8 MHz and the spectral resolution was 
62.5 kHz. The total bandwidth for the VLA D-array Hi 
observations was 6 MHz and the spectral resolution was 97.7 
kHz. Continuum subtraction weis done by spccifyiug line- 
free channels and using the MPS task UVLIN. The GMRT 
data were cleaned using the self-calibrated gains from the 
continuum data analysis. For the VLA data analysis, the 
bright continuum source M87 was subtracted using UVSUB, 



before the multi-channel data was continuum subtracted and 
CLEANed using IMAGR. 



3 RADIO CONTINUUM EMISSION 

3.1 A mini double- lobed radio source 

The VLA A-array image at 1490 MHz and the VLA B-array 
image at 4860 MHz with an angular resolution of 1.5 arc- 
sec (Fig. 2: left and middle panels) show the well-known 
doublc-lobcd structure of the radio source in the nuclear 
region (Hummel & Saikia 1991). The doublc-lobcd source 
has a total extent of ^12 arcsec (960 pc), with the western 
lobe being separated from the nucleus (the latter discussed 
in Section 3.2) by about 2.9 arcsec (230 pc) compared with 
8.9 arcsec (730 pc) for the eastern lobe. The spectral index 
image made from these scaled-array VLA images by con- 
sidering only those pixels which axe 5 times above the rms 
noise is also shown in Fig. 2 (right panel). The spectral in- 
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dex, defined as Soc v", varies from —1.5 to +0.5, while the 
error in the spectral index is ctq, — 0.08. The mean value of 
spectral index for the western lobe, the central region and 
the eastern lobe are —0.78, —0.84 and —0.61 respectively, 
the spectral index being steepest between the lobes. 

Given the steepness of a, the emission is clearly domi- 
nated strongly by synchrotron emission. However, it is im- 
portant to quantify what, if amy, contribution a thermal com- 
ponent might make. The latter would flatten the spectral in- 
dex from a steeper vahic, ont, to the flatter value, a, that is 
observed. In particular, we ask what the thermal contribu- 
tion would have to be in order to flatten the spectral index 
by an amount that is significantly outside of our error bar, 
ffa, in the spectral index map. We use the measured fiux 
densities of the two maps shown in Fig. 2, cut off to the 
same 5a level as the spectral index map. These values are 
Si.sGHz = 86.3 mjy and S4.86GHZ = 32.5 mjy, giving an 
observed total spectral index of a = —0.825. If the true non- 
thermal spectral index is ont = —0.825 — 2(Tq = —0.985, 
then (following Lee et al. 2001) we find that the thermal flux 
density required to alter the spectral index by this much is 
Si.5GHz{th) = 8.2 mJy. This corresponds to a required av- 
erage electron density of He — 7.2 cm~^ over a region whose 
equivalent spherical radius is 700 pc. While there may be 
small regions within the beam in which the density is higher 
than this (e.g. Kenney & Yale 2002), it is very unlikely that 
the electron density is this high, on average, over the ob- 
served region. Indeed, Machacek et al. (2004) find a mean 
density of rie « 0.03 cm~^ for the hot gas component over 
a roughly equivalent volume. We therefore conclude that 
the thermal contribution to these maps are negligible and 
continue with the minimum energy calculations under the 
assumption that a « unt- 

3.1.1 The minimum energy parameters 

The VLA images, which have the same spatial resolution and 
similar signal-to-noise (S/N) ratios, along with the spectral 
index map (Fig. 2) can be used to compute the following 
minimum energy parameters (see Pacholczyk 1970; Duric 
1991): the cosmic ray energy density ucr, the cosmic ray 
electron diffusion length Ld, the magnetic field strength B, 
and the particle lifetime t. This calculation can be done 
on a pixel-by-pixel basis assuming a particular geometry, as 
previously described in Irwin & Saikia (2003). 

In the case of the double-lobed source, we use a ge- 
ometry in which the line-of-sight depth is taken to be the 
measured average width of the mini double-lobed source, 
i.e. 4.21" = 344 pc with unity filling factor. The lower- and 
higher- frequency cut-offs of the spectrum are vi = 10' and 
1/2=10^^ Hz, respectively, and we adopt two different val- 
ues for the ratio of heavy— particle to electron energy, k=40, 
and k=100. The results are shown in Fig. 3 (a through d). 
In Fig. 3e, we also show a map of power, P=U/t, where U is 
the cosmic ray energy density mcr integrated along a line of 
sight and t is the particle lifetime. The result is identically 
P=(l-|-k)L, where L is the observed luminosity at a point. 
Thus the map of P closely resembles the map of total flux 
density but docs not match it exactly because the compu- 
tation of L requires an integration over frequency which is 
dependent on spectral index, and the spectral index is differ- 
ent at different points in the map. The map of P represents 



the rate at which cosmic rays must be accelerated in order 
to maintain equilibrium. 

Averages over the maps shown in Fig. 3 are given in 
Table 3. We have found that changes in the choice of line 
of sight distance and upper frequency cutoff make relatively 
small changes in the results in comparison to the choice of 
k (e.g. Irwin & Saikia 2003) . For the large range in adopted 
value of k, the results, overall, are within a factor of w 1.5 of 
each other. Note also, that this choice affects only the abso- 
lute scale of the maps and not the point-to-point variations. 
Beck & Krause (2005) have recently proposed adopting a 
system whereby the particle number ratio, rather than the 
energy ratio, is used in such calculations. For the mean spec- 
tral index of the mini double-lobed source (a ~ —0.8) the 
magnetic field values using their formalism are within a fac- 
tor of ~2 of ours (using the classical approach). As indicated 
in the table, the total power in the cosmic ray component 
is a few times lO'*'^ ergs s~^ and the total energy in cosmic 
rays is 10'''^' ergs. 

The estimates of the magnetic field, which determine 
the radiative lifetime of the particles (see below), are ~70 
/xG for k=40. It is of interest to compare this value with 
other galaxies which have radio bubbles or lobes powered by 
an AGN. Two of the well-studied galaxies with such features 
are NC6764 which has a total linear extent of 2.6 kpc (Hota 
& Saikia 2006) and Circinus which has an extent of 8.1 kpc 
(Elmouttie et al. 1998). The equipartition fields in these two 
galaxies for k=40 are ~25 and 45 /iG respectively. The value 
for NGC 4438 which is more compact is slightly higher. 

Of particular interest is the lifetime map (Fig. 3d) since 
it helps to determine whether or not the particles need to be 
accelerated in situ. We use lobe locations identified by the 
well-defined peaks in the power map (Fig. 3o), giving loca- 
tions for the western lobe of RA = 12*^ 27™- 45. "52 Dec -|-13° 
00' 32."9, and for the eastern lobe of RA = 12'' 27"" 46.n8 
Dec +13° 00' 26. "1. These positions correspond to projected 
separations from the nucleus of 2.6" (214 pc) and 9.2" (758 
pc) for the western and eastern lobes respectively. At these 
positions, the particle lifetimes are, for k = 40, t = 0.84 Myr 
(west) and t = 1.3 Myr (east) (somewhat lower if k = 100). 
Ignoring a possible linc-of-sight component to the velocity, 
the required particle velocity is only 250 km s^^ and 570 
km s~^ for the west and east lobes, respectively, if the rel- 
ativistic particles are supplied by the radio nucleus. These 
are lower limits, given the possible line-of-sight component 
to the velocity, but they arc sufficiently lower than c that 
in situ acceleration is not required, based on lifetime argu- 
ments. It is worth noting that the spectral index is typically 
flatter, the magnetic field stronger, and the particle lifetimes 
shorter at locations near the peaks of the two lobes in com- 
parison to regions between the nucleus and lobes. This be- 
haviour is typical of classical double-lobed sources where the 
peaks of emission are identified with regions where energy 
from the AGN is deposited. 



3.1.2 The eastern and western lobes/shells 

The higher resolution VLA A-array images at 4860 and 8460 
MHz, with angular resolutions of -^0.38 and 0.19 arcsec, 
respectively, are shown in Fig. 4. These reveal the promi- 
nent shell-like structure in the western lobe, noted earlier 
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Table 3. Minimum energy parameters." 



Model'' 


(eV cm~^) 


(pc) 


(mG) 


tf 
(Myr) 


UCR^ 

(10^5 ergs) 


Pcr" 
{W^^ ergs s-i) 


k=40 
k=100 


463 
775 


178 
137 


68.6 
88.8 


2.46 
1.67 


8.75 
14.6 


1.72 
4.24 



a See Pacholczyk (1970), Duric (1991) or Irwin & Saikia (2003). 

b A line of sight distance of 344 pc is adopted, k is the ratio of heavy particle to electron energy. 

c Cosmic ray energy density, averaged over the source. The total cosmic ray plus magnetic field energy 

density is 7/4 times these values. 
d Average diffusion length, 
e Average magnetic field strength. 
/ Average particle lifetime. 

g Total cosmic ray energy integrated over the source volume. 

h Total power, PcR = (1 +k) L, where L is the luminosity radiated by the electron component (see text). 



by Hummel & Saikia (1991), and also regions of emission 
from the eastern lobe. 

The structure of the western lobe, which is seen more 
clearly here than in the image of Hummel & Saikia, shows 
two main ridges of emission on the northern and southern 
sides with the peak of emission being on the northern ridge 
rather than the edge of the bubble farthest from the nucleus 
(Section 3.2). The western shell is also seen clearly in Ha and 
x-ray wavelengths (Kenney & Yale 2002; Machacek et al. 
2004). The eastern lobe is also known to exhibit a shell- like 
structure which is seen clearly at Ha and x-ray wavelengths 
(Kenney & Yale 2002; Machacek et al. 2004) although, in our 
high-resolution radio images, emission is seen only from the 
south-eastern and southern parts of the shell. The eastern 
shell, with a radius of ~0.5 arcsec (40pc), is smaller than 
the western one whose radius is 1.0 arcsec (80pc). 

Although there is a close correspondence between the 
shells seen at radio and Ha and x-ray wavelengths there 
are also significant differences. The eastern shell appears 
relatively more complete at Ha and at x-ray wavelengths 
whereas, as already noted, we have detected radio emission 
only from parts of the shell. Moreover, in the western lobe, 
both Ha and x-ray emission peak close to the nucleus which 
is rather weak at radio wavelengths. 
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3.2 The radio nucleus 

VLA A-array images at 4860 and 8460 MHz (Fig. 4) show 
clearly a compact radio source. The position of this source 
is RA 12'' 27"" 45. "67, Dec +13° 00' 31."54 at 4860 MHz, 
which agrees well with the position of RA 12'' 27"^^ 45.''66, 
Dec +13° 00' 31. "53 determined from the 8460-MHz image. 
These positions also agree with those of the optical and in- 
frared nuclei (Clements 1983; Keel & Wehrle 1993; Falco et 
al. 1999). The J2000 position of the infrared peak, which 
is likely to be least affected by extinction, is RA 12'' 27"* 
45.^67, Dec +13° 00' 31."54 (Keel & Wehrle 1993) and this 
value is also consistent with the measurements at optical 
wavelengths. 

The peak flux densities of the compact radio source esti- 
mated from our two images are 0.29 mjy/beam at 4860 MHz 
and 0.43 mJy/beam at 8460 MHz which yields an inverted 
radio spectrum with a spectral index of +0.69. The spectrum 
could be even more inverted if there is any contamination 
of the flux density at 4860 MHz by any extended emission. 



Figure 5. The GMRT image of NGG 4438 and 617 MHz with an 
angular resolution of ~33 arcsec. Contours: 4x(— 4, —2.82, 2.82, 
4, 5.65, 8 . . . ) mJy/beam. 

The highly inverted spectrum of this compact component 
suggests it to be the nucleus of NGC 4438 (see Section 5.1). 
It would be interesting to re-observe it for any possible signs 
of variability. The spectrum in this region appears steep in 
the low-resolution images (Fig. 2) because it is dominated 
by the extended bridge emission. The 8460- MHz image (Fig. 
4, right panel) shows evidence of a jet-like extension from 
the nucleus along a PA of ^^233° which is almost orthogo- 
nal to the orientation of the mini-double. The nature of this 
feature is unclear; a more sensitive image to trace its extent 
would be useful. 

3.3 The western extended radio emission 

In Fig. 5, we show the extended 617 MHz continuum emis- 
sion from NGC 4438, observed with the GMRT, over a fleld 
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Table 4. Radio flux density 



Frequency 


Stot. 


^double 


Reference 


MHz 


mJy 


mJy 




(1) 


(2) 


(3) 


(4) 


610 


324±30 




VTW2004; G2003 


617 


311±30 


151" 


p 


1400 


150±10 


86'' 


KE1983 


1420 


149±15 




VTW2004; G2003 


1490 




87 


P; HS1991 


4850 


97±9 




NKW1995; VTW2004 


4850 


70±10 




BWE1991 


4860 




33 


P; HS1991 


8600 


49±4 




VTW2004 


10550 


44±4 




NKW1995; VTW2004 



" Estimated from the peak flux density in our GMRT image. 
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References. BWE1991: Becker et al. 1991; G2003: Gavazzi ct al. 
2003; HS1991: Hummel & Saikia 1991; KE1983: Kotanyi & Ekers 
1983; NKW1995: Niklas et al. 1995; P: Present paper; VTW2004: 
VoUmer et al. 2004a. 

of view ~400 times larger in area than that shown in Fig. 2. 

The brightest emission, seen as an unresolved component at 
the centre of the map, contains the nucleus and the mini- 
lobe discussed earlier (Fig. 2). The western extended radio 
emission was initially reported by Allen et al. (1973) and 
Kotanyi, van Gorkom & Ekers (1983). The peak of the ex- 
tended emission is about --^1 arcmin from the nuclear region 
and extends over 2 arcmin from the nucleus of NGC 4438, 
or at least 10 kpc towards the south-west. 

The rms noise in Fig. 5 is 4 mjy/beam, which is higher 
than desirable, due to the presence of the strong radio source 
M87. This is also possibly responsible for the significant vari- 
ations in the total flux density values of NGC 4438 which 
have been reported in the literature (see Table 4). VoUmer, 
Thierbach & Wielebinski (2004a) have tried to determine 
systematically the radio continuum spectra of galaxies in 
the Virgo cluster and quote a spectral index of —0.67 for 
the total emission from NGC 4438 using a number of mea- 
surements between 600 and 10550 MHz. 

Subtracting the flux densities of the nuclear double- 
lobed source for which we have reliable measurements of the 
total flux density (Table 4) from the measurements of total 
flux density we can estimate the spectral index of the west- 
ern extended emission. Using the measurements of Vollmer 
et al. for the total flux density gives flux densities of 62 and 
64 mJy for the western extended emission at 1400 and ^5000 
MHz respectively. This yields a flat spectral index of '^O for 
the extended emission, even flatter than the value of ~— 0.5 
reported by Kotanyi et al. (1983). Adopting the value of 70 
mJy for the total flux density at 5000 MHz (Becker et al. 
1991) gives us a flux density of 37 mJy and a spectral index 
of ~— 0.4 for the extended emission. The western extended 
radio emission is visible in the 10550 MHz image of Niklas, 
Klein & Wielebinski (1995) possibly due to its flat spectrum. 
Estimating the flux density of the nuclear emission at 10550 
MHz by extrapolating it from lower frequencies and sub- 
tracting it from the total flux density also yields a similar 
flat spectral index of ~— 0.4 between 1400 and 10550 MHz. 
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© 
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Figure 6. The continuum flux densities of the total emission 
(filled circles), the inner double (open circles) and the western 
extended emission (filled triangles) of NGC4438 using the flux 
densities listed in Table 4. 

At low-frequencies the total flux density at 610 MHz 
(Vollmer et al. 2004a) is consistent with our estimate of 
311±30 mJy at 617 MHz with an angular resolution of 
46"x24" arcsec^ along a PA of 149°. Identifying the peak 
flux density of 151 m,Jy in the GMRT image at 617 MHz 
with the nuclear double gives a flux density of ~160 mJy for 
the extended emission. The peak flux density of 151 mJy is 
consistent with the extrapolated flux density of ~179 mJy 
at 617 MHz for the nuclear double using our VLA mea- 
surements. This yields a spectral index of ^—1.1 for the 
extended emission between 600 and 1400 MHz significantly 
steeper than the high-frequency spectrum (see Fig. 6) . This 
suggests that the low-frequency spectral index is dominated 
by non-thermal emission, while at high-frequencies contri- 
butions from thermal emission become important. This is 
a demonstration of both thermal and non-thermal radio- 
emitting material exisiting in the extra-planar gas. It is 
relevant to note that a flat spectral index of ~— 0.5 for 
extended extra-planar radio emission may also be due to 
re-acceleration of particles (e.g. Bell 1978; Duric 1986) as 
has been suggested for the Virgo cluster galax:y NGC 4522 
(Vollmer et al. 2004b). However, the difference in spectral 
index discussed above suggests that this is not the explana^ 
tion for NGC 4438. 



4 HI EMISSION 

Hi observations of NGC 4438 have been reported earlier with 
an angular resolution of 23x118 arcscc^ along PA=0° using 
the Westerbork telescope by Kotanyi (1981), and more re- 
cently with an angular resolution of ~40 arcsec (3.3 kpc) us- 
ing the VLA-D array by Cayattc et al. (1994) and again with 
the VLA D-array by Li & van Gorkom (2001). We observed 
this source with the GMRT with the goal of determining the 
structure of Hi in emission with higher resolution, and also to 
detect Hi in absorption towards the central compact source 
to identify any kinematic effects of the rmclcar activity on 
the Hi gas. Our highest resolution (see Table 5) was slightly 
higher than obtained in previous observations. To supple- 
ment our GMRT observations, we have also analysed VLA 
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Table 5. Observation log and observed parameters of the Hi-images. 



Tel. 


Obs. 


t 


Phase 


Seal 


Vcl 


Beam size 




map 


spec. 


Hi 




date 


min 


calib. 


Jy 


res. 


maj . 


min. 


PA 


rms 


rms 


flux 












km/s 


u 


// 


O 


m,Ty/b 


mjy 


Jy km/s 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


GMRT 


2002Aprl4 


300 


1254+116 


0.79 


13.2 


38 


35 


178 


1.2 


5.8 


2.9 


VLA-D 


1988Jul02 


330 


1252+119 


0.92 


20.7 


58 


48 





0.9 


2.0 


4.0 














127 


124 


84 


1.2 


4.6 


6.5 




-0.03 I ' ' ' ' ' ' ' 1 

-600 -400 -200 200 400 600 800 1000 

Heliocentric velocity (km/s) 

Figure 7. The global profile of Hi emission from the GMRT ob- 
servations with an angular resolution of ~36 arcsec (continuous 
line) and VLA D-array observations tapered to an angular reso- 
lution of ~125 arcsec (dashed line). 

D-array archival data to detect and study the disk emission, 
which has not been seen earlier, with high surface brightness 
sensitivity. We have also tapered the data to create another 
data set which emphasizes broadscale structure in order to 
detect and understand the properties of more extended Hi 
gas. The phase centre of this archival data set was the cen- 
tre of NGC 4438, unlike the observations reported by Li & 
van Gorkom (2001). The VLA results presented here axe 
of higher surface brightness sensitivity than those reported 
earlier and show several new features which we discuss in 
this paper. 

4.1 GMRT results 

We present here the global profile (Fig. 7) and the total- 
intensity Hi image from the GMRT observations (Fig. 8). 
These observations clearly show an elongated structure 
along a PA of ~35° which is very similar to the PA of 29° 
for the central stellar disk measured from molecular line ob- 
servations (Kenney et al. 1995). This elongated structure 
extends for ~6.5 kpc along the major axis in our GMRT 
image. It is roughly parallel to the stellar disk and ~4.1 kpc 
away from its mid plane. Earlier Hi images (e.g. Kotanyi et 
al. 1983; Cayette et al. 1994; Li & van Gorkom 2001) show 
the displacement of Hi gas on the western side of the stel- 
lar disk but these observations show the elongated structure 
more clearly. In addition to this elongated structure there 
are weak blobs of Hi emission towards the south. 

The Hi mass of the elongated feature is 1.5xlO^M0, 



^A¥g443B GMRT hi + D^lfiua^ 



13 a] a 




Figure 8. The moment zero image from the GMRT observations 
with an angular resolution of ~36 arcsec superimposed on the 
DSS blue-band image. Contour levels are (11.28, 16, 22.60, 32 
. . .)xlO^^ atoms cm~^, in steps of \/2. 

while the total H i mass estimated from the global pro- 
file (Fig. 7) obtained from the GMRT observation is 
~2x10^Mq. The global profile exhibits a sharp drop in in- 
tensity at ~— 40 km s^^ and a tail of emission on the red- 
shifted side. Although most of the H I gas appears red-shifted 
with respect to the systemic velocity of the galaxy (71 km 
s~^), the blobs towards the south appear blue shifted. The 
velocity gradient of the elongated structure shows that the 
north-eastern side is receding with a maximum heliocentric 
radial velocity of 107 km s~^ while the south-western side 
is approaching with a minimum radial velocity of —28 km 
s~^. These velocities are consistent with the Westerbork ob- 
servations reported by Kotanyi (1981). The middle point of 
the radial velocity range of the elongated feature from our 
observations is ~40 km s~^. The maximum line width in 
this elongated Hi emission structure is ~40 km s~^. The ob- 
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Figure 9. Left panel: The moment zero image from the VLA D-array observations with an angular resolution of ~53 arcsec superimposed 
on the DSS blue-band image. Contour levels: (3.84, 5.45, 7.69 . . .) x 10^^ atoms cm ^ in steps of V2. Middle panel: The corresponding 
moment one image showing the velocity field. Contour levels from south to north are -40, -20, 0, 20, 40, 60, 80 ... km s"! in steps of 
20 km s~^. Right panel: The corresponding moment two map showing the velocity dispersion. Contour levels: 20, 50 and 100 km . 



served sense of rotation of this feature is similar to that of 
CO emission seen in the stellar disk of the galaxy (Combes 
et al. 1988; Kenney et al. 1995; Vollmer et al. 2005), and 

lies close to the position of the extra-planar CO (1-0) omis- 
sion along a PA~20° (Vollmer et al. 2005). The extent of 
the extra-planar CO emission is very similar to that of the 
Hi feature seen here, although the peak of the CO emission 
lies slightly to the north-west while that of the Hl-elongated 
structure towards the centre of the feature. In this region the 
DSS blue band image shows the presence of dust plumes 
which appear very similar in nature to those seen in the 
Virgo cluster galaxy NGC 4402 (Crowl et al. 2005). 



4.2 VLA results 

We first present the results of the VLA D-array full- 
resolution observations with an angular resolution of 58x48 
arcsec'^ along a PA of 0°. Moment maps were made with a 
cutoff at 4(7 and integrated in velocity from —133 to veloc- 
ity -1-340 km s~^ (Fig. 9). With a higher surface brightness 
sensitivity than the GMRT observations, the total-intensity 
contours of Hi omission show that the elongated feature ex- 
tends further south with a total size of ~2' (9.8 kpc) along 
a similar orientation to that of the elongated feature seen 
in the GMRT observations. There is also Hi emission from 
the optical disk of the galaxy with the emission having a 
sharper edge towards the north and extending further to- 
wards the south with a similar total extent ~2' (9.8 kpc). 
There is a clump of Hi emission at the southern end of the 
disk emission with velocities in the range of ~— 20 to -|-20 
km s~ . The rnornerit map also shows a clump of Hi emis- 
sion towards the north between NGC 4438 and NGC 4435 
with heliocentric velocities ranging from -|-176 to +217 km 
s-i. 



The velocity field of NGC 4438 in this VLA D-array 
map is complex. In the central region of the disk within ~10 
arcsec of the nucleus the isovelocity contours are orthogonal 

to the central stellar disk with velocities ranging from 40 
km s~^ on the southern side to 240 km s~^ on the northern 
side, yielding a gradient of 10 km arcsec"^. At larger 
distances from the nucleus the isovelocity contours along 
the disk of the galaxy range from —80 near the southern 
clump to 280 km s~^ on the northern edge giving a veloc- 
ity gradient of 3 km s~^ arcsec"^. There is a hint of the 
isovelocity contours exhibiting a C-shaped curve or its re- 
flection towards the north-east and south-west regions of the 
extra-planar gas, which requires confirmation from observa- 
tions of higher spatial resolution. If confirmed, this would be 
reminescent of 'backward question marks' shaped isoveloc- 
ity contours in NGC 4522 (Kenney, van Gorkom & Vollmer 
2004) which is at a similar orientation and undergoing inter- 
action with the ICM of the Virgo cluster. The extra-planar 
gas in NGC 4438 has a smoother velocity field than the disk 
ranging from -|-10 on the southern side to 130 km s~^ on the 
northern side giving a lower velocity gradient of 1 km s~^ 
arcsec^^ over a similar length scale. It is worth noting that 
the blue-shifted velocities of —28 km s~^ in the southern 
side of the elongated feature seen in the GMRT and West- 
erbork observations is not apparent in the D-array image, 
suggesting that observations with different resolutions pick 
up sub structures with different velocities. 

Approximately 15 arcsec south of the nucleus centred at 
RA -12'' 27"^ 45.M, Dec 13° 00' 20", there is a small region 
of 15x10 arcsec^ along a PA~130° where the iso- velocity 
contours are closed and have values ranging from 40 to 80 
km s~^. These contours represent g£is which is approaching 
us relative to the gas in its vicinity. The moment 2 map 
shows that the velocity dispersion in this region is high with 
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NGC4438, the HI tail & IC3355 VLA-D 
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Figure 10. Hi-emission channel maps made with a resolution of ~125 arcsec. The contours levels are 1.2 x (—4, —2.82, 2.82, 4, 5.65 
. . .) mjy/beam. The systemic velocity of NGC 4438 is 71 km s~^. The velocity separation between adjacent channels is 20.7 km . 
The crosses in all the images denote the position of the radio nucleus of NGC 4438. The irregular galaxy, IC 3355, is located at RA 12'' 
26"" 51." 13, Dec +13° 10' 32."6 (see Sect.4.2.1). 



values in the range of 90—115 km s^^. The moment 2 map 
also shows that towards the north of this region, the velocity 
dispersion is typically 50—80 km s~^ over a total extent of 
^3.7 kpc approximately orthogonal to the PA of the disk. 
For the rest of the galaxy and the extra-planar emission the 
line widths are in the range of 20—50 km s~^ with the line 
widths being larger on the northern side. 

The VLA D-array data were also imaged with a IkA 
taper (resolution of ~ 125 arcsec) which shows more Hi gas 
primarily on the western side. The global profile (Fig. 7) 



superimposed on the one from the GMRT higher resolution 
observations shows that there is indeed more Hi gas towards 
the red-shifted side. The VLA spectrum also shows a sharp 
cut off on the blueward side of the spectrum as seen in the 
GMRT global profile. The total Hi flux density is 6.5 Jy km 
s~^ which is consistent with that of Arecibo mesurements 
by Giovanardi, Krumm & Salpeter (1983) who quote a line 
flux density of 6.1±0.6 Jy km s~^. Our VLA mesurements 
correspond to a total Hi mass of 4.4x10* M© which is ~2 
times larger than that estimated from the GMRT spectrum. 
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Figure 11. Left panel: The moment zero image from the VLA D-array observations tapered to an angular resolution of ~125 arcsec 
superimposed on the DSS blue-band image. Contour levels: (2.26, 3.22, 4.54, . . .)xlO^^ atoms cm^^ in steps of \/2. Middle panel: The 
corresponding moment one image showing the velocity field. Contour levels from south to north are —5, 20, 40, 60, 80, 100, 120, 140, 
160, 180 and 190 km s~^. Right panel: The corresponding moment two map showing the velocity dispersion. Contour levels are 20, 55 
and 70 km s"'^. 



The channel maps obtained with a spatial resolution 
of ~ 125 arcsec are shown in Fig 10. In addition to the Hi 
emission associated with NGC 4438 whose nucleus is marked 
with a cross in the figure, Hi emission is also seen at the 
position of the irregular galaxy, IC 3355 (see figure caption), 
and a large cloud of Hi, which we will refer to as the 'tail', is 
detected which is most prominent in the velocity channel at 
~— 9.5 km s^^. These features will be discussed more fully 
in the next two sections. 

Moment maps of NGC 4438 were made of the tapered 
VLA data with a cutoff at 4(t and integrated in velocity 
from -71 to 238 km s"^ (Fig. 11). The Hi total- intensity 
image with an angular resolution of 127x124 arcsec^ looks 
resolved with an extension towards the north. The northern 
clump seen in the full-resolution D-array image is coincident 
with this extension. In the tapered image the disk emission, 
the extra-planar elongated structure and the southern and 
northern clumps appear blended with more diffuse Hi gas. 

The velocity contours in the southern side of the disk arc 
approximately parallel to each other along a PA of ~— 80° 
while those on the northern side are oriented along a PA 
of ~— 45°. The velocity in the disk ranges from -|-20 to 140 
km s~^ with a velocity gradient of 0.9 km s~^ arcsec"^. 
The extra-planar gas, including the northern extension, has 
velocities ranging from ~20 to ~190 km s^^ yielding an 
average velocity gradient of ~0.5 km s^^ arcsec"^. At this 
resolution the Hi gas in the disk as well as the extended 
emission on the western side are rotating about a heliocen- 
tric velocity of -^110 km s~^ about an axis which is at a PA of 
—48°. There is some hint that the isovelocity contours of the 
extra-planar gas seen in this tapered image appear to curve 
inwards about this axis at larger distances from the disk, 
while in the full-resolution image they may diverge closer to 
the interface between the disk and the extra-planar gas, rem- 
inescent of the 'backward question mark' in the Virgo cluster 
galaxy, NGC 4522, which has been affected by ram-pressure 
stripping (Kenney, van Gorkom & VoUmer 2004). Some ev- 



idence of similar diverging iso- velocity contours can be seen 
in the galaxy NGC 2820, which is affected by ram pressure 
stripping due to the intra-group medium (Kantharia et al. 
2005) . It may be relevant to study the velocity structures in 
the wakes of ram-pressure stripped gas with features simi- 
lar to a von-Karman vortex street (e.g. Fig. 8 of Roediger, 
Briiggen & Hoeft 2006) to understand the above-mentioned 
iso-velocity contours. The moment 2 map shows that the 
width of the Hi line changes from ~20 km s~^ in the outer 
regions to ~87 km s~^ near the centre of the galaxy. 

It is also of interest to compare the Hi velocity field with 
those obtained at other wavelengths, although one must bear 
in mind that the Hi observations are of much coarser reso- 
lution than at other wavebands. The velocities of CO and 
Ha emission along PAs of 29° and 27° respectively, and Ha 
gas along a PA of 29°, all passing through the nucleus of 
the galaxy, have been presented by Kenney et al. (1995) 
and Chemin et al. (2005) respectively. The optical observa- 
tions having a seeing of 1.5—2 arcsec, the CO observations 
having an angular resolution of ~6 arcsec while the Hi im- 
ages have angular resolutions of ~50 and 120 arcsec for the 
full-resolution and tapered images respectively. As noted by 
Kenney et al. the CO and Ha velocities are symmetric and 
in good agreement within a galactocentric radius of ~10 arc- 
sec. Beyond this distance, the velocities become asymmetric 
with the CO velocities on the nortli-eastern side increasing 
to ~240 km s~^ at ~20 arcsec from the nucleus while the 
Ha velocities increase to ~275 km s~^ at ~40 arcsec from 
the nucleus. On the south-western side the CO extends to 
~15 arcsec with a minimum blue-shifted velocity of —105 
km s^^ but the Ha gas seen by Kenney et al. extends up 
to ~130 arcsec with a nearly constant velocity of ~— 20 km 
s~^. The velocities obtained by Chemin et al. (2005) are very 
similar to those of Kenney et al. on the north-eastern side, 
but goes to significantly higher blue-shifts of ^—150 km s~^ 
approximately 15 arcsec south-west of the nucleus. The Hi 
velocity field with a full-resolution of ~53 arcsec (Fig. 9) ex- 
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hibits a similar pattern with velocities ranging from 50 
km on the south-western side to ~260 km s~^ on the 
north-eastern side ~70 arcsec from the nucleus. 



4.2.1 IC 3355 

IC 3355 (VV 511, DDO 124, VCC 0945) is an irregular 
galaxy (Fig. 12) located -82 kpc from NGC4438. This 
galaxy is listed in the atlas of interacting galaxies by 
Vorontsov-Vel'Yaminov (1977) who also note the presence 
in its optical tail of '3 diffuse companions in a blue haze'. 
The B-band image shows a number of compact star-forming 
regions towards the eastern side of the galaxy. Spectroscopic 
observations at both radio and optical wavelengths listed 
in NED, the HyperLeda (http://leda.univ-lyonl.fr; Paturel 
et al. 2003) and GOLDmine (http://goldmine.mib.infn.it; 
Gavazzi et al. 2003) data bases and available information in 
the literature show the presence of two systems, one with a 
heliocentric velocity of ~— 10 km s~^ and the other at a he- 
liocentric velocity of ~162 km s~^ (de Vaucouleurs, de Vau- 
couleurs & Nieto 1979; Sulentic 1980; Haynes & Giovanelli 
1986; Hoffman et al. 1987; Cayatte et al. 1990; Huchra, 
Geller & Corwin 1995; Gavazzi et al. 2004). At radio wave- 
lengths Haynes & Giovanelli (1986) report the detection of 
a weak feature at 162 km s"'^ and note that the 'spectrum 
in severely blended with local Hi emission'. Our spectrum 
does not show any feature at —162 km s^^ (Fig. 12), and this 
velocity is perhaps incorrect. Hoffman et al. (1987) and Cay- 
atte et al. (1990) have reported the detection of Hi emission 
at a heliocentric velocity of —9 and —12 km s~^ respectively. 

The channel maps (Fig. 10) from the VLA D-array data 
clearly show the detection of Hi emission at the position of 
this galaxy, as noted in Sect. 4.2. The HI emission is centered 
at a position of RA 12*^ 26"^ 51. "1, Dec 4-13° 10' 33", which 
is coincident with the position of the optical galaxy: RA 
12*^ 26™ 51. "13, Dec +13° 10' 32. "6 (Yasuda, Ohamura & 
Fukugita 1995). In Fig. 10, Hi emission is detected in the 
channels with velocities of 11.1, —9.5 and —30.1 km s~^, 
but no emission is detected in the channels with velocities 
of 176.1 and 155.5 km s"'^, which are closest to the velocity 
system at 162 km s~^. 

At this resolution of ~125 arcsec the primary beam cor- 
rected total Hi flux density is 5.47 Jy km s~^ estimated from 
the global profile (Fig. 12, upper panel) and corresponds to 
a total mass of 3.72x10® M© for IC 3355. Note that this is 
almost as high as the Hi mass of NGC 4438, itself (4.4x10* 
Mq, Sect. 4.2), further emphasizing that the latter galaxy 
has been heavily stripped of its gas. A full-resolution VLA 
D-array image of IC 3355 (Fig 12, lower panel) shows the 
emission to be elongated approximately in the north-south 
direction, similar to that of the galaxy and consistent with 
the structure reported by Cayatte et al.(1990) with an angu- 
lar resolution of 21 x 17 arcsec^ along a PA of 75° obtained 
by combining VLA C— and D— array data. The western side 
of the Hi gas a sharp gradient but with the emission on the 
northern and southern extremeties bending marginally to- 
wards the west. The image presented by Li & van Gorkom 
(2001) detects extended emission on the western side which 
is not seen in either our image or that of Cayatte et al., pos- 
sibly due to lower sensitivity. The detailed structure of the 
Hi gas is possibly due to a combination of ram pressure and 
tidal interactions since IC 3355 along with NGC 4438, are 
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Figure 12. Upper panel: The primary beam corrected Hi emis- 
sion spectrum of IC 3355 made from a tapered VLA D-array im- 
age with an angular resolution of ~125 arcsec. Lower panel: The 
total-intensity Hi emission contour of IC 3355, observed with the 
VLA D-array with an angular resolution of ~53 arcsec superim- 
posed on the DSS blue-band image. Contour levels: (3.84, 5.45, 
7.69 . . .)xlOi^ atoms cm"^ in steps of V2. 



both likely to be part of the M86 sub-cluster, which appears 
to be merging with the M87 part of the cluster (see Kotanyi 
& Ekers 1983; Elmegreen et al. 2000; Schindler et al. 1999). 



4.2.2 An Hi-tail 

The channel maps at —9.5 and —30.1 km s^^ show extended 
diffuse emission which is seen most prominently at —9.5 km 
s~^ in Fig. 10. The global profile of the Hi tail is shown in 
Fig. 13 while the moment image showing the Hi emission 
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extending for ~10 arcmin is presented in Fig. 14. It is of in- 
terest to note that the Hi emission from the galaxy IC 3355 is 
also strongest in the velocity channel of —9.5 km s~^ which 
is close to its optical velocity of —10 km s~^, suggesting 
that the tail of Hi emission may bo of cxtragalactic ori- 
gin, although the possibility of it being Galactic foreground 
emission cannot be ruled out. If extragalactic, the tail is 
unlikely to be associated with the other companion galaxy, 
NGC 4435, which has a radial velocity of 801 km s"'^. In fact, 
the latter galaxy, although close in projection to NGC 4438, 
is likely physically farther away from NGC 4438, given the 
large (AV =730 km s~^) velocity difference between these 
two galaxies. 

Deep optical images of the Virgo cluster (Phillips & 
Malin 1982; Katsiyarmis ot al. 1998) show a faint stellar tail 
extending towards the north of NGC 4438. Deeper observa- 
tions trace the optical tail to a surface brightness of /Xv ~28 
mag arcsec"^ where it bends abruptly by ~90° to the west 
(Mihos et al. 2005). The 'knee' structure of the optical tail 
is expected in close and slow encounters in a cluster (cf. Mi- 
hos et al. 2005) and is unlikely to be caused by the interac- 
tion between NGC 4438 and 4435, and could be significantly 
older than ~100 Myr. The superposition of our Hi image on 
the deep optical image of Mihos et al. is shown in the bot- 
tom panel of Fig. 14. The Hi tail we have imaged is close in 
position and orientation to a significant part of the optical 
tail pointing northwards with the Hi extending till the bend 
in the optical tail. The northern tip of the tail is close to the 
edge of the half-power point of the primary beam; it would 
be useful to image the entire field to determine whether the 
Hi tail follows the optical one beyond the bend. This could 
also provide insights towards understanding whether the tail 
might be of cxtragalactic origin. 

Extragalactic Hi tails could be debris of tidal interac- 
tion, ram pressure stripped gas or intergalactic Hi in the 
cluster medium (sec Bekki, Koribalski & Kilborn 2005). The 
tail of Hi emission in NGC 4388 extends for ~110— 125 kpc 
and has a dynamical time scale of ~100 Myr (Oosterloo & 
van Gorkom 2005). The fine flux of the NGC4438 Hi tail in 
the velocity channel of —9.5 km s~^ is 2 Jy km s~^ yield- 
ing a mass of ~1.36xl0* Mq if it is at the distance of the 
Virgo cluser, which is about 30 per cent of the Hi mass of 
NGC 4438 or 36 per cent of the Hi mass of IC 3355. For 
comparison the mass of the Hi tail in NGC4388 is 3.4x10* 
M© (Oosterloo & van Gorkom 2005). 



5 DISCUSSION 

5.1 The radio-continuum emission 

Our detection of an inverted-spectrum nuclear component 
suggests that the small double-lobed radio source with a 
scale size of ~1 kpc is due to an AGN rather than a star- 
burst, consistent with earlier suggestions from x-ray observa- 
tions (Machacek et al. 2004). Although optically thick ther- 
mal emission from compact sources with sizes of ~5— 10 pc 
and electron densities of a few thousand cm~^ can have an 
inverted spectrum at cm wavelenghts, as seen in the compact 
star clusters in He 2-10 (Kobulnicky & Johnson 1999), the 
star formation rate in the nucleus of NGC4438 is modest 
(~0.05— 0.1 M© yr~^) and estimates of the electron den- 
sity in the nuclear region are small, ^10 cm~^ (e.g. Kenney 
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Figure 13. The primary beam corrected Hi emission spectrum 
of the Hi tail made from a tapered VLA D-array image with an 
angular resolution of ~125 arcsec. 



et al. 1995; Machacek et al. 2004). Prom a compilation of 
the structures of nearby galaxies with a starburst and/or 
an AGN, Hota & Saikia (2006) have suggested that such 
bubble- or lobe-like radio structures are more likely to be 
seen in sources with an AGN rather than a starburst. 

One of the most striking aspects of the structure on 
scales of ~1 kpc is the radio lobes and that the western 
lobe is clearly seen as a shell-like structure. The shell-like 
structures are also seen in Ha and x-ray wavelengths on 
opposite sides of the nucleus and are closely related to the 
radio structure. The Ha and x-ray emission possibly arise 
in regions where the bubbles of synchrotron-emitting plasma 
interact with the surrounding ISM. 

The radio lobes are very asymmetrically located with 
the ratio of separations of the eastern lobe from the rmclcus 
to that of the western one being -^3, while the correspond- 
ing flux density ratio is ~0.60 and 0.57 at 1.4 and 5 GHz 
respectively. A higher density on the western side where the 
ISM of NGC 4438 appears displaced could provide a viable 
explanation of the asymmetry in location relative to the nu- 
cleus (cf. Kenney & Yale 2002). The lobe on the denser side 
is also expected to be more luminous due to a higher ef- 
ficiency of conversion of beam energy into radio emission, 
and better confinement of the radio lobe (cf. Eilek & Shore 
1989; Jcyakumar et al. 2005). However, this simple scenario 
does not explain why the western lobe, which is closer to the 
nucleus and brighter, is also larger in size by a factor of ~2. 

A possible solution to this inconsistency comes from 
the results of numerical simulations of propagation of jets. 
Hydrodynamic simulations of light, large-scale jets in a de- 
creasing density profile show that the jet bow shock under- 
goes two phases, first a nearly spherical one and secondly 
the well-known cigar-shaped one (Krauso 2002; Krause & 
Camcnzind 2002; Carvalho & O'Dea 2002). The shell-like 
structure of the western lobe is suggestive of the first phase 
of the development of the bow-shock. In this scenario, the 
western jet has not yet entered the cigar phase, and deposits 
its radio-emitting plasma in a bigger part of the bubble, al- 
most filling the region within the bounds of the bow shock. 
On the other hand, the eastern jet appears to be in the cigar 
phase and should therefore have a fairly regular backflow 
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Figure 14. The total-intensity contours of Hi emission of the pos- 
sible tail and the galaxies NGC 4438 and IC 3355 superimposed 
on the DSS blue-band image (upper panel) and the deep optical 
image (lower panel) from Mihos et al. (2005). The Hi contours 
levels are (1.13, 1.60, 2.27 . . .)xlOi'' atoms cm-^ in steps of V2. 



around it, which flows back into the central parts diffusing 
and mixing with the shocked external gas. A similar expla- 
nation has been suggested for the more distant radio galaxy 
with a starburst, 3C459, which exhibits a similar asymme- 
try in the location and structure of the lobes (Thomasson, 
Saikia & Muxlow 2003). 

For a galaxy moving with a velocity of ~1000 km s~^ 
through the Virgo cluster, ram pressure could also affect 
the observed structure of the source. Using plausible values 
from Combes ct al. (1988) and Vollmcr ct al. (2005) the ram 
pressure is approximately few times 10~^^ dynes cm~^. For 
comparison the pressure in the eastern and western lobes 
arc X^IO^^" dynes cm^^, which corresponds to a cosmic ray 
energy density of ~200 eV cm^'^ (sec Fig. 3a). This is signif- 
icantly larger than the ram pressure, while the pressure in 
the interlobe region is comparable to that of ram pressure. 
The lobes of emission are unlikely to be significantly affected 
by ram pressure due to the motion of the galaxy. 

The extended emission displaced to the western side 
of the disk of the galaxy has a spectral index of — 1.1 
between 600 and 1400 MHz significantly steeper than the 
high-frequency spectrum which is likely to be between 
and ^—0.4. This suggests that the low-frequency spectral 
index is dominated by non-thermal emission, while at high- 
frequencies contributions from thermal emission become im- 
portant. This extended emission is unlikely to be due to 
the AGN. It is also interesting to note that radio emission 
has not been detected from the disk of the galaxy; the ex- 
tended emission possibly represents plasma which has been 
displaced from the disk along with other components of the 
ISM. 



5.2 The Hi gas 

The GMRT and VLA observations show the structure of 
Hi emission from both the disk and extra-planar gas asso- 
ciated with NGC 4438. Hi emission has also been detected 
from IC 3355 as well as a 'tail' of Hi gas extending for ~10' 
between NGC 4438 and IC 3355. 

The extra-planar Hi gas in NGC 4438 has an approx- 
imately linear structure with a total size of ^9.8 kpc and 
displaced from the disk by ~4.1 kpc. The elongated Hi struc- 
ture overlaps with the CO extra-planar emission but has a 
larger length than the CO feature which has been repro- 
duced in simulations involving both tidal interactions and 
ram pressure stripping by the ICM of the Virgo cluster with 
ISM-ISM interactions between NGC 4438 and NGC 4435 
playing a relatively minor role (VoUmer et al. 2005). There 
is a hint that the isovelocity contours in the extra-planar 
gas appear to initally diverge while at larger heights they 
curve inwards. These features, which require confirmation 
from higher spatial resolution observations with adequate 
surface-brightness sensitivity, are reminescent of the galaxy 
NGC 4522 also in the Virgo cluster, which has been inter- 
preted to be due to ram-pressure stripping. The VLA D- 
array full-resolution image shows a region of closed velocity 
contours where the gas appears to deviate from the overall 
rotation and approaches us. A similar feature is also seen in 
the Ha observations at a similar location but with a much 
higher approaching velocity (Chemin et al. 2005). There is 
an x-ray source at this location and it is tempting to spec- 
ulate that the outflowing gas may be related to energy re- 
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leased from the x-ray source. The high velocity dispersion 
of the Hi gas in this location would be consistent with such 
an interpretation. 

The faint Hi-tail we detect towards the north of 
NGC 4438 is close to a stellar tail seen earlier in deep opti- 
cal observations. The velocity of the Hi tail is 10 km s~^ 
similar to that of Hi emission from IC 3355. Although these 
aspects suggest that the tail may be in the Virgo cluster, the 
possibility that it may be a Galactic feature cannot be ruled 
out. If it is at a distance similar to that of NGC 4438, the 
Hi-tail has a total extent of ~50 kpc and a mass of 1.4x10* 
Mq . It is relevant to note that in NGC 4388 which has an Hi 
tail which extends for ~110 kpc (Oostcrloo & van Gorkom 
2005), there is no evidence of optical light along the length 
of the tail except for faint streamers close to the disk. Al- 
though such Hi tails could be debris of tidal interaction, ram 
pressure stripped gas or intergalactic Hi in the Virgo cluster 
medium (see Bekki, Koribalski & Kilborn 2005), the pres- 
ence of optical emission light would suggest that it is likely 
to be due to tidal interactions. It would be interesting to 
make deep Hi observations of NGC 4438 to image this tail 
beyond the 'knee' of the optical tail and also find counter- 
parts of the other faint optical tails and arc-like structures 
listed by Mihos et al. (2005). 



galaxy for the first time as well as the detailed velocity struc- 
ture of the disk and extraplanar gas. The cxtraplanar gas 
appears displaced to the west of NGC 4438 by ~4.1 kpc and 
has a mass of 1.8 x 10* Mq which is ~1.5 times higher than 
the mass in the disk. The extraplanar gas appears to bo in 
regular rotation and may be rotating more slowly than the 
gas in the disk. The extraplanar gas appears to be rotating 
about a velocity of ~110 km s~^, which is redshifted com- 
pared with the systemic velocity of the galaxy. The velocity 
structure is possibly affected by ram pressure and is rcmines- 
cent of the classic ram pressure stripped galgixy NGC4522 
in the Virgo cluster. 

Our observations also reveal a large tail of Hi extending 
north of NGC 4438. Its velocity is very similar to that of 
the irregular companion galaxy IC 3355. Also the Hi tail 
is close to that of a faint optical tail seen in deep images, 
suggesting that the tail may be in the Virgo cluster although 
the possibility that it might be Galactic emission cannot be 
ruled out. If it is located close to NGC 4438, its mass is 
1.36 X 10* Mq, it has a total size of ~50 kpc and extends 
for ~75 kpc from NGC 4438. The possible association of 
starlight suggests that it may have formed tidally due to 
interactions with one of the galaxies in the M86 merging 
sub-cluster. 



6 CONCLUSIONS 

We have made the most extensive radio continuum and Hi 
observations of the Virgo Cluster galaxy, NGC 4438, and its 
environment over a variety of spatial resolutions using the 
GMRT and the VLA. We have presented new data as well as 
archival data which have been imaged using improved data- 
reduction parameters, revealing new features related to this 
galaxy. The radio continuum data span spatial resolutions 
of 0.2 to 33 arcsec and the Hi data span spatial resolutions 
from 35 to 125 arcsec. 

From our high-resolution, radio-continuum data, it is 
clear that the galaxy harbours an AGN whose position 
agrees with the optical and IR nucleus and which has an 
inverted radio spectrum. This confirms ealier speculations, 
based on x-ray data, that an AGN is present and is consis- 
tent with the presence of radio lobes in the galaxy which 
are more likely to be associated with an AGN rather than a 
starburst (Hota & Saikia 2006). 

We have presented the clearest view of the miru-double- 
lobed radio source in NGC 4438 and presented the first maps 
of the minimum energy parameters in the radio lobes. Note 
that this is one of the few spiral galaxies known that contains 
a clear double-lobed radio source. The western lobe extends 
to ~ 230 pc from the nucleus and the eastern lobe to ~ 730 
pc from the nucleus. The lobes breaJj up into shells at higher 
resolution, visible in x-ray and Ha observations as well. In 
our radio data, the western shell is clearly delineated. 

At lower resolution, a large radio continuum extension 
is visible, from the 617 MHz GMRT data, on the south- 
western side of the galaxy. This feature extends as far as 
10 kpc from the nucleus of the galaxy. An analysis of the 
spectral index of this feature suggests that in addition to the 
largely non-thermal emission it may contain a contribution 
from thermal gas. 

The Hi observations show emission from the disk of the 
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